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^ i Abstract 

In a microwave background polarization map that covers only part of the sky, it is 
impossible to separate the E and B components perfectly. This difficulty in general 
makes it more difficult to detect the B component in a data set. Any polarization map 
' can be separated in a unique way into "pure E," "pure B" and "ambiguous" compo- 

nents. Power that resides in the pure E {B) component is guaranteed to be produced 
by E (B) modes, but there is no way to tell whether the ambiguous component comes 
from E or B modes. A polarization map can be separated into the three components 
either by finding an orthonormal basis for each component, or directly in real space by 
, using Green functions or other methods. 
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Detailed characterization of the polarization of the cosmic microwave background (CMB) ra- 
diation will be one of the main frontiers in cosmology in the coming years. The first detection 
of CMB polarization |2| and the first measurement of the large-angle correlation between 
polarization and temperature E] have already taken place, and numerous experiments in 
the present and near future El 13 El E] will surely provide a wealth of polarization data. 
^ ■ A CMB polarization map consists of maps of the Stokes parameters Q and U. 1 The Stokes 

parameters are of course coordinate-dependent objects: Q+iU is a spin-2 field. A key insight 
in CMB polarization theory is that the most natural way to express such a spin-2 field is to 
decompose it into scalar and pseudoscalar pieces, generally called E and B [TDl[TTJ . This E/B 
decomposition is crucial in analyzing polarization data. In particular, scalar perturbations 
(such as density variations) produce only E-type polarization (in linear theory), so detection 
of a B component could provide evidence for vector or tensor perturbations. The search 
for the B component is therefore likely to be of enormous importance: this component is 
capable of telling us about tensor (gravity-wave) perturbations produced during inflation 
[T2"| 113 j . probing the inflationary epoch far more directly than any other observations. 

In standard models, the B component is considerably weaker than the E component, and 
so is likely to be difficult to detect [Hj. To make matters worse, the E/B decomposition is 
unique only for a full-sky map. This means that in the absence of complete sky coverage, 

1 No circular polarization is expected, so we ignore the Stokes parameter V. 
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Figure 1: Left panel: An E Fourier mode. Center panel: A B Fourier mode. Right panel: 
An E "hot spot." 

unless one is very careful, E/B confusion can reduce the detectability of the B component 

damn. 

One way to quantify the problem of E/B leakage is to observe that the space of all possible 
polarization maps over any given region of sky can be decomposed into three orthogonal 
subspaces: a "pure E" subspace, in which all power is guaranteed to come from E modes, 
a "pure B" subspace, in which all power is guaranteed to come from B modes, and an 
"ambiguous" subspace, in which there is no way to tell whether the power came form E or 
B. This decomposition was worked out in spherical harmonic space for a spherical cap in |17j . 
and the general formalism is described in [TH] . In the latter work, explicit recipes are given 
for finding orthonormal bases ( "normal modes" ) for all three subspaces as eigenfunctions of 
the bilaplacian operator on the observed region. Furthermore, in the case of a pixelized map, 
we give an efficient way of finding approximately pure and ambiguous modes by solving a 
discrete eigenvalue problem. 

An understanding of the _E/-B/ambiguous (hereinafter E/B/A) decomposition is likely 
to be of great importance in designing future experiments. For instance, the presence of am- 
biguous modes significantly increases the optimal sky coverage for a degree-scale 5-detection 
experiment [To] . 

2 E and B modes 

We begin by summarizing some fundamental properties of E- and 5-type polarization maps. 
For further details, see |18j . Suppose for the moment that we have a polarization map that 
covers a small enough patch of sky to allow the use of the flat-sky approximation. In this 
case, it is natural to work in Fourier space, where the separation into E and B is quite 
simple. As shown in Figure ^ a Fourier mode is an E mode if the polarization direction 
is parallel / perpendicular to the wavevector and is a B mode if the polarization direction 
makes a 45° angle with the wavevector. Any E polarization map is a superposition of such 
E Fourier modes; an example is the E "hot spot" shown in the right panel of Figure ^ (To 
get a B spot, simply rotate the polarization by 45° at each point.) Realizations of E- and 
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Figure 2: Gaussian random maps. One is an E map; one is a B map. For the reader who 
wishes to test his or her ability to tell E from B, we reveal which is which in the references 

EH]- 



.B-type Gaussian random fields are shown in Figure 121 

The Fourier-space description of the E/B decomposition is very helpful in developing an 
intuitive understanding of the mixing of E and B modes in an incomplete or pixelized map. 
In a map that covers only part of the sky (say of size L), we can achieve resolution ~ L^ 1 
in k space. That is, our estimate of a Fourier mode with wavevector k will actually include 
contributions with a range of wavevectors centered on k with a spread of ~ L . Since 
wavenumbers with different directions will be mixed together, and since the direction of the 
wavevector is crucial to the E/B separation, mixing is inevitable, with the worst problems 
occurring on the largest scales probed (k ~ L~ x ). 

In a pixelized map, problems crop up on the small-scale end as well, due to aliased 
power. When a mode with a wavevector k beyond the Nyquist frequency is aliased to a 
lower frequency, it typically is mapped into a mode whose wavevector points in a completely 
different direction. As a result, aliased power has nearly complete E/B mixing. For this 
reason, avoiding aliasing (by oversampling the beam) is even more important in polarization 
experiments than in temperature experiments. 

In real space, the two types of polarization maps satisfy differential equations on Q and 
U. Any E map must satisfy 

dL-p = o, (1) 

where P = is the polarization map, and the differential operator Dg is 2 

°- =(<£%)■ < 2 > 

2 These equations are written in the flat-sky approximation for simplicity. The general equations are given 
in [H]. 



"vector" ( ^ ) , while the conjugate operator T) B turns a polarization vector into a scalar. 



Note that the operator D E acts on a scalar function to produce a two-component polarization 
U 

Similarly, any B map must satisfy • P = 0, with 

The easiest way to see that these equations are true is to verify that they work for arbitrary 
E and B Fourier modes. 

In building intuition about E and B modes, it is extremely helpful to bear in mind the 
analogy with vector fields: an E mode is the spin-2 analogue of a curl-free vector field, and 
a B mode is the analogue of a divergence-free vector field. This means that the operators 
Db and T) E are like the curl and divergence respectively. 

This analogy immediately suggests a conjecture. Any curl-free (divergence-free) vector 
field can be written as the gradient (curl) of a potential. Perhaps a corresponding rule works 
for spin-2 fields. This conjecture turns out to be correct: any E map can be written as 

P E = B E ij E (4) 

for some scalar "potential" ip E . Similarly for B maps: P E = D^. The reason all of 
this works is that ■ Db = ■ D# = 0, the analogues of the familiar vector identities 
V-Vx = V x V = 0. We note for future reference the other useful identity 

T> E = B B D B = (V 2 ) 2 = V 4 , (5) 

the bilaplacian. 3 For vector fields instead of spin-2 fields, the laplacian shows up on the right 
instead of the bilaplacian: V • V = V x Vx =V 2 . 4 

Incidentally, just as in the case of vector fields, potentials are guaranteed to exist only 
when the functions are defined over a simply-connected region. If the observed patch of sky 
has "holes" in it, then there are E modes that cannot be derived from a potential. Figure El 
shows an example. Despite its strikingly B-\ike appearance, this is in fact an E mode (i.e., 
D E applied to it gives zero) , although it cannot be expressed as D E applied to a scalar 
potential. 5 The E/B/A decomposition described in jTHj still works on such regions: modes 
like this one show up automatically as ambiguous modes. 



3 If we do not make the flat-sky approximation, this becomes V 2 (V 2 + 2). 

4 The astute reader may have noted a slight problem with this equation: it doesn't appear to be true. In 
three dimensions, there is an extra term in the V x Vx equation. The extra term is not present, however, in 
the case of interest here, where the operator is being applied to a (pseudo)scalar function in two dimensions. 

5 At the risk of belaboring the obvious, the vector-field analogue of this is the magnetic field of a long 
straight current, v = eg/r in polar coordinates. This is curl-free over any region not containing the origin, 
but if the region completely surrounds the origin it cannot be expressed as the gradient of a potential. 
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Figure 3: An E mode that cannot be expressed as the gradient of a potential. 

3 E/B/A Separation 

We begin by summarizing some key results from ^H]. Any biharmonic function, i.e., any 
function if) such that 

VV = 0, (6) 

generates a pair of polarization maps D ^if> and D gifj that simultaneously satisfy the condi- 
tions for E and B modes. We call such modes "ambiguous." Furthermore, all ambiguous 
modes can be expressed in this way. 

We define a pure E mode to be one that is orthogonal to all B modes (including the 
ambiguous modes, which are after all B modes). Every pure E mode can be written in 
the form Dg?/), where if) satisfies both Dirichlet and Neumann boundary conditions on the 
boundary of the observed region. (That is, both if> and the normal component of vanish 
on the boundary.) These conditions imply that the pure E modes always have polarization 
parallel or perpendicular to the boundary at the edges of the map. 6 Similarly, any pure B 
mode can be written as D^. Pure B modes always hit the boundary of the map at a 45° 
angle. 

A natural way to generate an orthonormal basis of pure E and B modes is to find 
the eigenfunctions of the bilaplacian that satisfy both Dirichlet and Neumann boundary 
conditions. (The bilaplacian has a complete set of such eigenfunctions.) For the case of a 
disc, a sample of ambiguous, pure E, and pure B modes is shown in Figure 0] 

In some cases, it may not be convenient to find a complete basis of normal modes; it may 
be preferable to perform the E/B/A separation on a map directly in real space. We present 
here a brief sketch of some ways this can be done. 

First, suppose that we have a polarization map P that covers the entire sky (taken for 
simplicity to be a plane rather than a sphere). We can express P as the sum of an E piece 

6 Contrary to our statement in [18] . this is not a sufficient condition. In addition to being parallel / 
perpendicular on the boundary, there is an another constraint that must be satisfied for an E mode to be 
pure. 
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Figure 4: The top panel shows examples of ambiguous modes for a disc. The bottom panel 
shows examples of pure E modes. To generate pure B modes, rotate the polarizations in 
the pure E mode maps by 45° at each point. 




(7) 




where the E Green function must satisfy ■ G#(x) = 5(x) and ■ G# = 0. Explicit 
forms for the Green functions are easily found: 



in polar coordinates. These Green functions give the E and B "response" to any given point 
in a polarization map. For instance, if a map contains a delta-function spike in Q at the 
center, then the E and B maps will be as shown in Figure |SJ In this Figure and the following, 
because of the large range of polarization magnitudes, we use a logarithmic greyscale rather 
than the lengths of the lines to indicate magnitude. The polarization maps plotted in this 



figure are D^(— cos26'/47r) and D#(sin 29 /An), the results of substituting P(x) = 





(9) 




into equation (JHJ). 
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Figure 5: The E and B maps that result from decomposing an input map consisting of a 
delta function in Q. In this figure and the next, the lines, which are all of equal length, 
indicate direction. The magnitude of the polarization is indicated by the gray scale, which 
is logarithmic and covers a factor of 10 4 . 

This figure illustrates the nonlocal character of the E/B decomposition. Note that the 
Green functions for the potentials do not go to zero at large distances: to find the potentials, 
we need to specify P arbitrarily far away. Since the actual polarization is a second derivative 
of the potential, however, the response in the E and B polarization maps to a delta-function 
impulse does decline as the inverse square of the distance. 

In the case where the map covers only part of the sky, the same approach can be used to 
get the pure E and B components of the map. However, the Green functions must be replaced 
by functions that satisfy the appropriate Dirichlet and Neumann boundary conditions. Such 
functions can in principle be found, but in practice this is unlikely to be an efficient way to 
perform the E/B/A separation. 

A much more efficient way to perform the E/B/A separation is to split the task up into 
two steps: 

1. Separate the map into E and B components without worrying about purity. 

2. "Purify" each of the two maps by projecting out the ambiguous component. 

The first step can be done in several ways, the most efficient being to Fourier transform 
the maps and do the separation mode by mode. (Equation (|SJ) can also be used, but this 
will generally be slower.) If the map is of an irregular shape, it can be padded out to a 
convenient rectangular array in any way you like before Fourier transforming; the resulting 
E/B decomposition will be correct. Different ways of doing this padding will in general lead 
to different decompositions, but they will differ only in where the ambiguous modes show 
up. 

One way to perform the second step is to find yet more Green functions. Suppose we have 
an E polarization map that we wish to purify. We first find a potential ipE that generates 
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Figure 6: The E and B maps produced by the "ambiguous" Green functions g\ and g 2 for 
a disc. As in the previous figure, the gray scale indicates the magnitude of the polarization 
covers a factor of 3 x 10 5 logarithmically. 

this map. (This is also easily done in Fourier space.) We then subtract a biharmonic function 
from this potential so that the residual has the correct (Dirichlet and Neumann) boundary 
conditions. One way to do this is 

^| ure (x) = ip E (x) - f dip [V>i?(^)£i(x; ip) + {h- V)^0)£ 2 (x; <p)]. (10) 
Jan 

Here Q represents the observed patch of sky. We use ip to label points on the boundary dQ. 
The two Green functions g\ and g 2 are biharmonic functions (i.e., generators of ambiguous 
modes) that satisfy delta-function boundary conditions. Specifically, 

g 1 (ip' ] ip) = «% - ip'), (n ■ V)g x (<pf\ cp) = 0, (11) 

and the reverse for g 2 . For the case of a disc, these Green functions can be calculated 
analytically and are shown in Figure El 

The two two Green functions give the corrections that must be applied to an E polar- 
ization map to purify it of any failure to meet the correct boundary conditions at any given 
point. (Since there are two boundary conditions, in general two corrections must be applied). 
Both produce polarization maps that decrease like with increasing distance: one is inversely 
proportional to the distance cubed, while the other is inverse square. This confirms what has 
been noted elsewhere [T71 [TH] : the ambiguous modes tend to be largest near the boundary. 

In practice, equation(|l(Jj) is unlikely to be the most efficient way to purify a map. Given 
the potential ipE fo r an E map that we wish to purify, it will generally be much faster to use 
other numerical methods to find a biharmonic function %pA whose value and first derivative 
match ^ on the boundary. Subtracting T)e^a from the original E map will purify it. 
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4 Conclusions 



The presence of an ambiguous component in a CMB polarization map significantly affects 
the science that can be derived from the map. In particular, experimenters attempting 
to detect B modes in degree-scale experiments should be sure to take E/B confusion into 
account when designing experiments |16j . 

In a pixelized map, aliasing of small-scale power significantly worsens the problem of 
E/B mixing. Especially considering the extremely blue spectrum predicted for the E modes, 
experimenters searching for B modes should be sure to heavily oversample the beam. 

The E/B/A decomposition is likely to be useful in analyzing data sets. Strictly speaking, 
it is not necessary to perform any such decomposition: it is possible in principle to com- 
pute the likelihood function L{Cf , Cj 3 ) for the E and B power spectra directly from the 
raw Q and U maps. However, decomposing the map first may increase the efficiency of the 
analysis, especially if we are willing to simply throw away the ambiguous modes: the likeli- 
hood function will then factor, with the each pure subspace's contribution depending only 
on the corresponding power spectrum. Aside from any increase in efficiency in evaluating 
likelihoods, the E/B/A decomposition will be useful in checking for systematic errors and 
foreground contaminants and for purposes of visualization. 

In [T^j, we presented methods for finding orthogonal bases for the E/B/A components. 
In many cases it may be preferable to perform the decomposition directly in map space, 
without finding a basis. The Green function approaches presented here are unlikely to be 
numerically efficient for large data sets; their primary purpose is to aid in the visualization 
of the E/B/A decomposition. The decomposition can be performed fairly efficiently in map 
space, however, in the manner briefly sketched at the end of the last section: perform an 
"impure" E/B decomposition first, and then purify each component by finding a biharmonic 
function satisfying appropriate boundary conditions. Numerical methods for efficiently find- 
ing biharmonic functions exist. 

Finally, let us note that the E/B decomposition shows up in the analysis of weak lensing 
data {e.g., |2TH I21|). The methods described herein may be useful in such analyses. 
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